Spatially-resolved probing of a non-equilibrium superconductor 
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Spatially resolved relaxation of non-equilibrium quasiparticles in a superconductor at ultra-low temperatures 
was experimentally studied. It was found that the quasiparticle injection through a tunnel junction results in 
the modification of the shape of the I-V characteristic of a remote 'detector' junction. The effect depends on 
the temperature, injection current and proximity to the injector. The phenomena can be understood in terms 
of creation of quasiparticle charge and energy disequilibrium characterized by two different length scales Aq« 
~ 5 flm and Aj> ~ 40 flm. The findings are in good agreement with existing phenomenological models, while 
more elaborated microscopic theory is mandatory for detailed quantitative comparison with experiment. The 
results are of fundamental importance for understanding electron transport phenomena in various nanoelectronic 
circuits. 

PACS numbers: 74.25.Fy, 74.40.+k, 74.50.+r, 74.78.Na 



I. INTRODUCTION 

Conversion of electric current at an interface between dif- 
ferent materials is a common process in any electric circuit. 
At nanoscales the whole system might behave as an inter- 
face if the dimension(s) are comparable to the characteristic 
relaxation length. Of particular interest are boundaries with 
a superconductor where electric current converts from single 
electrons to Cooper pairs. In a superconductor at a finite tem- 
perature there are always non-paired electrons called equi- 
librium quasiparticles. In the presence of additional distur- 
bance their concentration can be increased by non-equilibrium 
quasiparticles. If these excitations originate solely from pair- 
breaking they equally populate excitation spectrum with mo- 
menta pq > pf (electron-like) and pq < pp (hole-like), where 
Pf is the Fermi momentum. Such a perturbation contributes 
to the longitudinal (or energy) mode (L) contrary to transverse 
(or charge imbalance) mode (T) where the two branches of the 
excitation spectrum are non-equally populated. 1 Within the 
Keldysh formalism one can define the non-equilibrium dis- 
tribution functions /l,t- 2 The usual electron distribution func- 
tion is 2f = 1-fL-fr- 

These phenomena have attracted attention in mid-70s re- 
sulting in an impressive number of papers. 3,4 Those early 
experiments were mainly performed on sandwich-type flat 
structures not adequate for spatially-resolved studies. Agree- 
ment between experiment and theory was established reliably 
mainly in the high temperature limit T — > T c . The opposite 
limit T -C T c has been poorly investigated. Recently the in- 
terest to the problem increased, 5 " 8 while the understanding 
is still far from being satisfactory. In addition to the general 
interest in physics of superconductivity, relaxation of quasi- 
particle excitations is a bottleneck process in many nanoelec- 
tronic applications: solid state coolers, 9 " 11 cold/hot electron 
bolometers, 12 and Cooper pair boxes. 13 Experimental study 
of the corresponding phenomena is the subject of the paper. 

The paper is organized as follows. In the next Section, 



we describe our setup and the experimental conditions. In 
Sec.|III]we discuss our experimental data. We also formulate 
a phenomenological model which is in a good agreement with 
our findings. The model enables us to obtain the energy and 
charge modes relaxation lengths. Finally, we summarize the 
results in Sec. [IV] 



II. EXPERIMENT 

Multi-terminal nanostructures (Fig. [TJ were fabricated us- 
ing electron beam lithography and UHV evaporation of alu- 
minum (S) and copper (N) separated by naturally-grown alu- 
minum oxide (I). The sample layout is conceptually close to 
the layout, described in Refs. |5| and [lj. Electrons were in- 
jected into the superconductor through a large-area NIS injec- 
tor with typical tunnel resistance ~ 3 kil. The induced dis- 
equilibrium was measured by remote detector junctions with 
typical tunnel resistance ~ 50 k£l. Even at the lowest values 
of the utilized tunnel resistances one can neglect the proximity 
effect. 

The 3 He 4 He dilution refrigerator was located inside the 
electromagnetically shielded room with only battery-powered 
analogue pre-amplifiers inside. Typical measuring set-up 
(Fig. [TJ top panel) consisted of two circuits: injector and de- 
tector. For detection two complementary configurations were 
used: either a single NIS junction, or two connected in se- 
ries NISIN junctions from both sides of the superconduc- 
tor. The latter was found to be less sensitive to dc offsets 
presumably originating from parasitic thermoelectric effects. 
Care was taken not to increase the temperature of the system 
by the Joule heating at the injector by making the electrode 
rather massive (Fig.QJ bottom panel). An extra NIS junction 
(marked 'Thermometer' in Fig. [TJ and located ~1 jim away 
from the injector junction) was used to measure the electron 
temperature of the injector rj(/,) by fitting the experimental 
h{VT,Ii = const) dependencies with the text-book expression 
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FIG. 1 : (Color online) Top: schematics of the sample and measure- 
ments. Electron temperature of the normal metal injector is measured 
with the NIS junction 'Thermometer'. The quasiparticle relaxation 
in the superconducting long bar is determined by the 'Detectors': ei- 
ther a single NIS junction, or a pair of connected in series NISIN 
junctions from the opposite sides of the superconductor. Bottom: 
scanning probe microscope image of the Si/SiO A - substrate with two 
closest detector pairs and the injector (left most massive probe). Typ- 
ical thickness of the aluminum wire, copper detectors, and injector 
are 25 , 40 and 80 nm, and the width 400, 180 and 1600 nm, respec- 
tively. Critical temperature of the aluminum wire is about 1.35 K and 
the low temperature mean free path t ~20 nm. 



for the tunnel current of a NIS junction. Here and below 7, (V,), 
Id(Vd), and It(Vt) are the I-V characteristics of the injector, 
detector, and thermometer tunnel junctions, respectively (see 
Fig-E top panel). Once the relation 7^(7,-, T — const) was es- 
tablished, the corresponding 'Thermometer' circuit was dis- 
connected for the rest of the experiments. 



III. RESULTS AND DISCUSSION 

To avoid confusion it is instructive to summarize the defi- 
nitions of various parameters of the dimensionality 'temper- 
ature' hereafter used in the paper. The T stands for the bath 
temperature of the system measured by the two RuO A resis- 
tors thermally anchored to the mixing chamber and the mas- 
sive copper sample holder. The RuO r sensors were calibrated 
against the nuclear orientation thermometer, and during the 
experiments their readings differ less than by few mK. 

The T e ' and T~ denote the electron temperatures of the nor- 
mal metal injector and detectors, respectively. The T e ' and 
Tg are determined by fitting the corresponding Uiyi) and 
Idtyd) experimental dependencies by the well-known equilib- 
rium expression for the tunnel current of a NIS junction. The 
T e ' depends on the injection current, 7^(7,), due to the triv- 
ial Joule heating of the normal metal injector. At the lowest 
bath temperature T ~20 mK and the highest injection currents 



Ii ~ 1 jiA the electron temperature of the injector increased by 
8Tg =T* — T~10Q mK. The electron-phonon interaction is 
rather weak at ultra-low temperatures, 9 and the superconduct- 
ing 'body' of the sample, thermally anchored to the substrate, 
has low thermal conductivity. Hence, the phonon tempera- 
ture of the remote detectors should not deviate significantly 
from the bath temperature T. On the contrary, the electron 
temperature Tf of the normal detectors, determined by fit- 
ting the zero-injection 7^(^,7, = 0) dependency, was always 
higher than the base temperature T. The effect originates 
from the poor thermal coupling typical for ultra-low temper- 
atures and from inevitable electron heating coming from the 
noisy EM environment. At a base temperature T ~ 20 mK 
and zero injection 7, = the absolute value of the detector 
electron temperature offset ST* = T~ — T varied from 10 to 
40 mK depending on the coupling of the particular junction 
with the measuring circuit. As the detectors are decoupled 
by the two NIS junctions from the remote 'hot' injector, it is 
reasonable to assume that for the given detector its electron 
temperature T* should not depend on the injection current /,■: 
T e d (I i >0,T)~Tf{I i = 0,T). 

The Tg denote the superconductor electron temperature. 
Formally, it should be determined taking into account a com- 
plicated energy balance. 9 From the very basic considerations 
one may expect that at non-zero quasiparticle injections 7,- > 
the actual (thermodynamic) temperature of the superconduc- 
tor Tg, which enters into the non-equilibrium distribution 
function fj(E,Tg), should be higher than the bath temper- 
ature T. However, within the phenomenological formalism 
employed in the paper, Tg cannot be determined. More elab- 
orated microscopic analysis (to our best knowledge, currently 
absent) should be used to determine the non-equilibrium (and 
- essentially asymmetric with respect to the chemical poten- 
tial) distribution function. Otherwise, any arbitrary value of 
Tg substituted into a symmetric (e.g. equilibrium) distribu- 
tion function fj(E, T e s ) gives exactly the same tunnel current 
of a NIS junction. 

Finally, the T* characterizes the reduction of the supercon- 
ducting gap A. When the superconducting gap is suppressed 
by the quasiparticle injection one can use two complemen- 
tary descriptions: either to deal with the bare experimen- 
tal data A (/,■), or to convert the non-equilibrium energy gap 
into an effective temperature T* using conventional (equilib- 
rium) BCS expression. In the latter case T* just means: what 
would be the equilibrium temperature of the superconductor 
T = T* resulting in the corresponding reduction of the gap 
A(r/,7, > 0) = A B cs{T*Ji = 0). It should be emphasized that 
T* is only a convenient parameter of dimensionality 'temper- 
ature' and has no direct relation to the (non-equilibrium) ther- 
modynamic temperature Tg. As it will be shown below, the 
impact of quasiparticle injection on a superconductor is man- 
ifold. In particular, the longitudinal (energy) mode cannot be 
described solely by the reduction of the gap, or alternatively - 
the effective temperature T* . 

Injection of non-equilibrium quasiparticles into a supercon- 
ductor results in a deviation of its density of states (DOS) Ns 
and distribution functions f[ T from equilibrium. Under cer- 
tain assumptions the distribution function can be found by de- 
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FIG. 2: (Color online) (a) Sample C8 Id(Vd) characteristics at various 
injections for the double junction NISIN detector located at L,j=0.8 
flm away from the injector. Inset: the corresponding dlj/dVj de- 
pendencies. Slight asymmetry of the dl c i /dVj at high injections pre- 
sumably originates from the not perfect equality of the each of the 
NIS junctions connetced in series, (b) Sample C2 NIS detector lo- 
cated at Lij=2. 8 flm experimental dl^/dVj dependencies atT=18±l 
mK at three injections eV,/A ; = (O), 10 (A) and 20 (□), and at 
T = 200 and T = 400 (0) mK at zero injections. Solid lines are 
the theoretical best fits with the fitting parameters indicated in the 
plot. 



convolution of the tunneling I- V dependencies. 15 In more gen- 
eral cases, self-consistent solution of Keldysh-Usadel equa- 
tions is required. 2 For the interpretation of the results we 
take the simplified approach postulating equilibrium func- 
tional forms of the DOS and distribution function fs, but with 
parameters deviating from their equilibrium values, 



N S (E, A,r) =| Re{ (E + iT/2)/d(E + iT/2f - A 2 } 



fs(E,T e s ^*) = 1/ {exp ((£ - ll*)/k B T e s ) + 1 } . 



(1) 

(2) 



Deviation of f[ from equilibrium formally corresponds to 
an increase of r/ above the bath temperature T, while de- 
viation of fj - to a non-zero shift of chemical potential fi* . 
Considering the above expression for fs with two parameters, 
Tg and jj.*, one might account for both the longitudinal and 
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FIG. 3: (Color online) Deparing parameter T vs. injection energy for 
three detectors at L (d =0.8 (A), 2.8 (□) and 7.8 (*) flm. Left inset: 
sample C2 NIS detector located at L id = 7.8 flm r(/,) vs. A(/,-) at 
various injection currents /, at T = 18 mK (O) and T = 200 mK 
(A). Right inset: sample C2 T vs. injection energy at L,y=2.8 flm at 
T =17.5 (□), 200 (■) and 400 mK (ffl). 



the transverse mode excitations. For the analysis of the ex- 
perimental data we utilized the simplified phenomenological 
approach with three fitting parameters: DOS broadening T, 16 
superconducting gap A, and effective quasiparticle chemical 
potential /x*. 17,18 All parameters can be defined by fitting ex- 
perimental Id(Vd) or/and dI c i/dV c i dependencies. It should be 
noted that the utilized phenomenological analysis, based on 
the postulated equilibrium functional form of the distribution 
function fs(E,Tf disables determination of the (thermo- 
dynamic) temperature of the superconductor 7jr , More elab- 
orated model (to our best knowledge - currently absent) is 
mandatory to deduce the non-equilibrium fs(E, Tf,ji*) and 
the corresponding 7 1 /. 

As the injected non-equilibrium quasiparticles relax on a 
certain length, the fitting parameters should depend on the in- 
jection rate (energy eVi or current /, = VjG' NN ) and the distance 
between the detector and the injector L^. The Id{V c i) depen- 
dency of a NIS detector in the presence of non-equilibrium 
injection is 



Id = (G d NN /e) / N S (E,AX) [fr - fr] dE, 



(3) 



where Gf/ N is the normal state conductivity of the junc- 
tion, J*t (E,V d , 7f ) = tanh((£ + eV d )/2k B T e N ) - tanh((£ - 
eVd)/2kBTf) the equilibrium value of fj in the normal de- 
tector, and fj,(E,Li,i,Tg) the local non-equilibrium value of 
fr m the superconductor. The first term corresponds to the 
text-book equilibrium expression for a NIS tunnel current, ex- 
plicitly depending on properties of the superconductor only 
through its DOS N S (E, A,F). 

The shape of the experimental I c i(V t i) characteristic depends 
on the quasiparticle injection (Fig. [2j». By fitting I c i(V c i) or/and 
did /dVd dependencies one finds that the parameters T and A 
depend on the injection rate, proximity to the injector L,y, and 
temperature T (Figs. [3j |4j. The smearing of the Ij(Vd) de- 
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FIG. 4: (Color online) (a) Superconducting gap A vs. injection en- 
ergy atL id =0.S (A), 2.8 (□), 7.8 (★) and 17.8 (0) flm. Inset: reduc- 
tion of the superconducting gap measured by the same NIS detector 
at L id =7.$ /im at T =17.5 (open star), 200 (*) and 400 (crossed 
star) mK. (b) Effective temperature T* vs. injection energy for sam- 
ples C3 (18 mK, open symbols), C2 (19 mK, solid symbols) and C8 
(24 mK, crossed symbols) at L id =Q.% (A), 2.8 (□) and 7.8 (*) [lm. 
Lines are guides to the eye. Inset: spatial decay of T* averaged over 
several samples at eV,/A,= 70 (□), 16 (O) an d 2.5 (A). Lines are 
fits assuming exponential dependence ~ exp(— L^/Aj*). 



pendence [rounding of the 'corners'at eV c ( ~ A (/,•)] can be as- 
signed to the parameter r(/,) which enters the superconduc- 
tor DOS and is physically related to the finite quasiparticle 
lifetime. 16 For the same injection rate /; the temperature de- 
pendence r(T, /, = const, Ljj = const) (right inset in Fig. [3} 
most probably originates from the more intensive relaxation 
of the quasiparticles (i.e. from the reduction of the quasipar- 
ticle lifetime responsible for finite F) due to inelastic scatter- 
ing on phonons. The total number of quasiparticles entering 
the superconductor per unit time is proportional to /,-, while 
only the fraction ~ exp(— L^j \/Dxr) reaches the detector. 
At low temperatures the recombination time Xr can be es- 
timated from the relation, To/tr ~ (r/r c ) 1//2 exp(— A/fc^r), 
with To being the electron-phonon scattering time. 19 ' 20 These 
simple arguments qualitatively explain the spatial, energy, 
and temperature dependencies of Y (Fig. |5J, while a more 
elaborated model, presumably incorporating the effect of the 



environment, is required for a quantitative analysis. 

Neglecting the superconducting phase, the gap equation for 
the pairing potential A, to be solved self-consistently with the 



Keldysh-Usadel equations, is A — X Jq f[{E 7 Li l t)Re(F)dE, 
where X is the electron-phonon coupling constant, F the 
anomalous pair amplitude, and f[(E,Ljd) the local non- 
equilibrium value of /l in the superconductor.— At any en- 
ergy E the function f[(E,Lid) is smaller than its equilibrium 
value tanh(E /2ksT), 22 reducing A (Fig. HH). The effect is 
smaller at the lowest temperatures and increases with increas- 
ing bath temperature (Fig.|4^, inset). Utilizing the well-known 
BCS temperature dependence for the superconducting gap one 
can assign an effective temperature T*, which would provide 
the same effect on the energy gap as the quasiparticle injec- 
tion: A{T e s ,I i ,L id )=A BC s{T*,Ii = 0,L id ) (Fig. |4p). Assum- 
ing an exponential decay r*=r*(0)exp(— L,y/Ar*) one finds 
the characteristic length A T * = 40 jxm ± 20 jXm (Fig.HJ?, in- 
set). The large uncertainty comes from the weak Agcs(T) 
dependence at low temperatures and the limited set of Ly 
points (number of detectors). Once again we would like to 
remind that T* is a convenient parameter introduced to ac- 
count for the reduction of the superconducting gap A(/,) and 
has no direct relation to the thermodynamic temperature of 
the superconductor Tg which should enter the corresponding 
non-equilibrium distribution function. It might even happen 
that in the limit of strong injection and weak quasiparticle- 
quasiparticle scattering Tg cannot be defined at all, while the 
non-equilibrium energy gap A (/,■), or alternatively T*(Ii), can 
be still deduced from the experimental Id(Vd) characteristics. 

We believe that the reduction of the gap A(/,) and the in- 
crease of r(/,) are the two manifestations of the same phe- 
nomenon - energy disbalance (longitudinal mode) accompa- 
nying injection of the non-equilibrium quasiparticles. Presum- 
ably, in a comprehensive microscopic description (to our best 
knowledge, currently absent) the longitudinal and the trans- 
verse modes may be entangled both contributing to the A(/, ) 
and r(/, ) dependencies. Within the 'orthodox' BCS model the 
(temperature-dependent) smearing of a NIS detector Id(Yd) 
characteristic comes solely from the distribution function of 
the normal electrode f" (E,Vd,Tf ) and not from the DOS 
or distribution function of the superconductor. Temperature 
of the superconductor Tf enters solely through the temper- 
ature dependence of the gap A(r/). In a typical equilib- 
rium situation T = Tg = Tg. Employed in the paper anal- 
ysis deviates from that 'orthodox' model by introduction of 
the Dynes smearing r[/,] entering the DOS of the supercon- 
ductor Ns(E, A [/,■], r [/,-]), Eq. ([T), at a given, but not defined, 
temperature of the superconductor Tg, which indeed leads to 
the extra (non-temperature dependent) smearing of the I c i(Vd) 
dependence, Eq. ([3]). 

The employed analysis of the experimental data based on 
the equilibrium expression for the tunnel current [first term 
in Eq. ©J with the DOS deviating from the 'orthodox' BCS 
model [Eqs. (fl}, (O] enables separation of the two impacts 
A(/, ) and r(/,-). The reduction of the gap A(/ ( ) (Fig. |4]l leads to 
a 'sharp' I c i(Vd) characteristic: as if the temperature is still low 
Tg ~ T, but the energy gap is reduced A (7, > 0) < A(7, = 0). 
The impact of the finite r(/,) manifests itself as the smear- 
ing of the Id(Vd) characteristics at the 'corners' eVd — ±A(/,) 
which cannot be accounted for any reasonable value of T® . 

It is important to emphasize that the experimental non- 
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FIG. 5: (Color online) Sample C3 NIS detection, (a) Zoom of the 
sub-gap region of the l d {Vd) characteristic at different injection volt- 
ages V;=0, 3, 5, 10 and 15 mV. Lines correspond to the same injec- 
tions of reverse polarity, (b) Excess current I e f vs. injection energy 
eVi/Ai at different distances L,y. Lines are fits utilizing Aq* = 3.6, 
3.6, 5.2 and 6.4 flm for the detectors located at L id = 2.8, 7.8, 17.8 
and 37.8 flm, respectively. Inset: spatial decay of the excess current 
measured at | eVj/A; |=33.7 (O), 6.7 (□) and 4.5(A). (c) Normal- 
ized chemical potential fl* vs. injection energy measured at various 
temperatures by the same detector at L,y=7.8 fim. 



equilibrium /^(V^,/,- >0,r) dependencies cannot be described 
by the 'orthodox' BCS model just assuming that the temper- 
ature of the system has increased above the bath tempera- 
ture: T T e N , T 7 / > T. At T < T c the shape of the exper- 
imental Id(Vd,Ii > 0,T) characteristics is qualitatively differ- 
ent from the reference equilibrium characteristics IdiVchU = 
0,T') taken at elevated temperatures T' > T (Fig. Wp)- Those 



zero-injection high-temperature Id(Yd) characteristics can be 
nicely fit by the equilibrium BCS expression for the tunnel 
current [first term in Eq. ([3])] assuming small intrinsic Dynes 
broadening of the DOS, Eq. ©, with r(/, = 0)/A(7, = 0) ~ 
0.02 (Fig. in a reasonable agreement with existing data 
on aluminum. 21 However, one runs into a contradiction trying 
to fit the experimental Id (Vd , /,■ > 0,T) characteristics using a 
single fitting parameter 7/,, and the first term in the expression 
(01 for the NIS tunnel current. To properly account for the 
gap reduction A(/,) one should set the best fit 'equilibrium' 
temperature 7},-, = T* 3> T (compare Figs. [2J} and[4j)) being 
significantly higher than the temperature of 'hottest' part of 
the system - the injector. While to account for the observed 
broadening at the gap edge, the effective fitting temperature 
Tfj t should be set significantly lower 7V& ^ K <C (Fig. Wp)- 
Summarizing, one should conclude that at non-zero injections 
the Id(Vd,Ii > 0,T) characteristics cannot be described by a 
single 'equilibrium' temperature The two independent 
parameters, A(/ ( ) and r (/,■), are necessary. The dependencies 
A(/,) and r(/ ; ) on the injection rate are different. At ultra- 
low temperatures T <C T c and modest injections eVj/Ao <10 
the suppression of the gap is rather small (Fig. |4^), while 
the increase of the DOS broadening parameter T is quite pro- 
nounced (Fig. [3]). At elevated temperatures the dependence of 
r on the injection rate is rather weak (insets in Fig. [3). We 
would like to emphasize that the necessity to introduce the two 
independent parameters A(/,) and r(/, ) is the result of the em- 
ployed phenomenological approach based on Eqs. (H])-©. We 
would wish to believe that in an elaborated microscopic model 
dealing with the essentially non-equilibrium distribution func- 
tion fs(E,T^) a single parameter - the truly thermodynamic 
temperature of the superconductor 7jr - would be sufficient to 
account for the whole variety of the experimental data. 

A remarkable feature of the tunnel current expression is 
the existence of the non-zero excess current — Id (Vd = 0) 
originating from the independence of fj(E,Ljd,T^) on the 
bias volatage Vd [second term in Eq. (01]. The effect has 
been observed in NISIS' flat sandwiche s 23 ! 24 and recently 
in multiterminal NIS structures. 5 Within our phenomenolog- 
ical approach, the experimentally measured excess current 
is linked to the effective quasiparticle chemical potential, 24 
If = (G d NN /e) /i*. A typical example of a NIS detector sub- 
gap Id{Vd) characteristic is shown in Fig. [5^. The excess 
current increases with increase of the injection rate revers- 
ing its sign with the change of the injection polarity. The 
effect decreases with increasing distance to the injector Lid 
(Fig. I5J3) and the bath temperature (Fig. [5};). Slight asym- 
metry of If with respect to zero injection (Fig. [5^) probably 
originates from a parasitic residual thermo-electric potential. 
It should be noted that in the NISIN detector geometry the 
excess current is absent since the contributions from the two 
NIS junctions connected in series and located at the same dis- 
tance Lid (from the opposite sides of the superconducting bar, 
see Fig. [U, cancel each other. The double NISIN junctions 
were used only as a complementary configuration in the ex- 
periments on the energy mode relaxation (Fig. [2^). 

The solution of the diffusion equation for the charge imbal- 
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ance Q* gives the excess current 

where N{0) = 1.08 x 10 47 1/Jxm 3 is the normal state DOS 
for aluminium at the Fermi level, and F* is a smooth function 
varying from at no injection to 1 at <?V, > > A,. 24 Substituting 
the mean free path £ ~ 20 nm, Fermi velocity \>f = 1 .36 x 10 6 
m/s, diffusivity D — 1 /3vf£ and cross section <7 = 25 nm x 
400 nm one gets an acceptable fit for the dependence of the 
excess current on the injection energy (Fig. |5t>). The best-fit 
charge imbalance length Ag* = ^/Dig* varies from 3.5 to 6.5 
/im being in a reasonable agreement with recent findings. 5 ' 8 

Weak temperature dependence of the charge imbalance re- 
laxation Q* ~ I? ~ jj,* (Fig. [5};) at ultra-low temperatures is 
expected: in the absence of an effective electron-phonon inter- 
action the remaining (weak) mechanism is the elastic scatter- 
ing due to the gap inhomogeneity, or/and anisotropy or even 
just a finite supercurrent. 24 Contrary to the transverse mode 
relaxation, the relaxation of the longitudinal (energy) mode al- 
ways requires the presence of inelastic mechanism(s). At the 
ultra-low temperatures of the experiment ^I<A there are 
very few equilibrium phonons capable to provide the quasi- 
particle energy relaxation. The only (potential) source of the 
phonons with energies above the bath temperature T is the 
Joule heating in the injector. However, how it has been already 
discussed, even in the worst-case scenario of the strongest 
injection the diffusion of heat from the 'hot' injector into 
the superconductor is not sufficient to account for our find- 
ings. The maximum rise of the electron temperature in the 
injector ~ 100 mK. Even assuming a perfect electron- 
phonon coupling in the injector, the energy of the correspond- 
ing phonons is still much smaller than the energy > 2A re- 
quired for a pair of hole-like and electron-like quasiparticles 
to form a Cooper pair. However, as the quasiparticle recom- 
bination finally does happen on 'astronomically' large (for a 
superconductor) scales At* — 40 fim, some weak mechanism 
should exist. For example, emission / re-absorption of non- 
equilibrium phonons or photons, might provide the required 
energy relaxation. Certainly, these mechanisms should also 
contribute to the relaxation of the transverse mode. The ob- 
servation Ag* < At* supports our conclusion that at ultra-low 
temperatures the elastic scattering is an important channel of 
the charge imbalance relaxation. 

Substantial quantitative difference between At* and Ag* at 
ultra-low temperatures, T T c , is in a qualitative agreement 
with previous calculations . 19 ' 20 The explicit condition when 
one can neglect the elastic scattering has been obtained in 



Ref.lH (T c - T)/T c < (T £ r c .)- 2/3 , where T E is the energy 
relaxation time (measured in 1/K). Analysis of this expres- 
sion leads to the conclusion that the only physical limit when 
always Ag* ^ Ag is at T — > T c , where Ag* diverges. Formally 
it means that in the normal state there is no branch relaxation 
mechanism. Our experiments were made at very low temper- 
atures much smaller then T c . At such low temperatures, to our 
best knowledge, there is no general criterion to a priori relate 
Ag* with At* . 

IV. CONCLUSION 

In conclusion, we performed spatially resolved measure- 
ments of the non-equilibrium quasiparticle relaxation in su- 
perconducting aluminum at ultra-low temperatures T <C T c . 
The analysis employs the phenomenological interpretation 
postulating equilibrium functional dependencies of DOS and 
distribution function of the superconductor, while assuming 
that the broadening parameter T, effective chemical potential 
jj.* , and superconducting gap A depend on the rate of quasi- 
particle injection and the proximity to the injector. Two scales 
A T * = 40 ±20 fim and Ag* = 5 ± 1.5 fim are responsible 
for relaxation of the energy (longitudinal) and charge (trans- 
verse) modes, respectively. Both quantities were measured 
simultaneously in the same samples using the same experi- 
mental technique eliminating sample / measurement artifacts. 
This is the central result of the paper. It should be emphasized 
that both the energy and the charge disequilibrium are univer- 
sal phenomena which should be taken into consideration in 
a broad class of systems involving electron, 9 " 13 spin and/or 
coherent non-local 26 transport. Despite the reasonable agree- 
ment with the phenomenological description, application of 
the relaxation approximation for the essentially spatially inho- 
mogeneous problem is not fully justified at low temperatures 
r<T ( , A more elaborated microscopic model is required for 
a quantitative analysis. We hope that our findings will trigger 
the corresponding research activity. 
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